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Saturn’s moon Enceladus harbours a global’ ice-covered water ocean”*. The Cassini 
spacecraft investigated the composition of the ocean by analysis of material ejected 
into space by the moon’s cryovolcanic plume* ®. The analysis of salt-rich ice grains by 
Cassini’s Cosmic Dust Analyzer” enabled inference of major solutes in the ocean 
water (Na+, K*, Cl, HCO, , CO,” ) and its alkaline pH°™. Phosphorus, the least abundant 
of the bio-essential elements” “, has not yet been detected in an ocean beyond Earth. 
Earlier geochemical modelling studies suggest that phosphate might be scarce in the 
ocean of Enceladus and other icy ocean worlds’. However, more recent modelling of 
mineral solubilities in Enceladus’s ocean indicates that phosphate could be relatively 
abundant”. Here we present Cassini’s Cosmic Dust Analyzer mass spectra of ice grains 
emitted by Enceladus that show the presence of sodium phosphates. Our observational 


results, together with laboratory analogue experiments, suggest that phosphorus 

is readily available in Enceladus’s ocean in the form of orthophosphates, with 
phosphorus concentrations at least 100-fold higher in the moon’s plume-forming 
ocean waters than in Earth’s oceans. Furthermore, geochemical experiments and 
modelling demonstrate that such high phosphate abundances could be achieved in 
Enceladus and possibly in other icy ocean worlds beyond the primordial CO, 
snowline, either at the cold seafloor or in hydrothermal environments with moderate 
temperatures. In both cases the main driver is probably the higher solubility of 
calcium phosphate minerals compared with calcium carbonate in moderately alkaline 
solutions rich in carbonate or bicarbonate ions. 


Enceladus’s global ocean lies under an ice crust and above a rocky core 
where tidal dissipation’*”’ is suspected to drive hydrothermal activ- 
ity*®. There are several lines of evidence describing how volatile’ and 
dissolved materials?*” from the rocky core are either emitted by the 
plume in the gas phase or incorporated into ice particles, respectively. 

Cassini’s Cosmic Dust Analyzer (CDA) recorded time-of-flight (ToF) 
mass spectra with a mass resolution of m/Am = 10-50 (refs. 10,20) for 
cations generated by high-velocity impacts of individual grains onto 
the instrument’s rhodium target. The E-ring of Saturn is formed by ice 
grains escaping Enceladus’s plume into orbits around Saturn”, and 
hence the analysis of these grains by CDA provides important insights 
into the composition of the subsurface ocean—including a rich variety 
of organic compounds® °—with much better statistics compared with 
data from the rare occasions when Cassini traversed the plume itself. 

In a previous analysis it was shown that a fraction of these ice 
grains, called Type 3 particles, contain salts at significantly higher 
concentrations than all other ice grains in the plume and the E-ring*”. 
By co-addition of spectra from 107 individual detections inthe E-ring, 


average salt concentrations of 0.5-2.0% by weight (around 0.07- 
0.30 M) have been deduced for Type 3 grains, with NaCl, NaHCO,, 
Na,CO, and KCl as the most abundant constituents. Because it was 
inferred that these grains are derived from aerosolized ocean water’, 
the salt chemistry of the plume and E-ring grains would reflect that 
of the ocean. 

Phosphorus is an element essential for planetary habitability 
but to date it has not been detected in an ocean beyond Earth. Previous 
geochemical modelling suggested that phosphate might be scarce in 
the ocean of Enceladus and other icy ocean worlds®"*. However, more 
recent modelling of mineral solubilities in Enceladus’s ocean suggests 
that phosphate could be relatively abundant (roughly 1077-107 M), 
depending on the compositional characteristics of the ocean such as 
its pH and carbonate content, as well as on attainment of certain chemi- 
cal equilibria between ocean water and seafloor alteration minerals”. 

In this work we present CDA mass spectra of a population of E-ring 
ice grains that show the presence of sodium phosphates. We then per- 
form laboratory analogue experiments to quantitatively establish 
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Fig. 1| CDA cation spectrum co-added from nine baseline-corrected 
individual ice grain spectra. The mass lines signifying a high-salinity Type 3 
spectrum are Na’ (23 u) and (NaOH)Na’ (63 u) with secondary Na-rich signatures 
of (H,O)Na* (41 u) and Na,* (46 u). Sodium phosphates are represented by 
phosphate-bearing Na-cluster cations, with (Na,;PO,)Na* (187 u) possessing 

the highest amplitude in each spectrum followed by (Na,HPO,)Na’ (165 u) and 
(NaPO,)Na* (125 u). The first two unlabelled peaks at the beginning of the 
spectrum are H* and C*, stemming from target contamination”. a.u., arbitrary 
units. 


that Enceladus’s ocean is rich in dissolved phosphate. Lastly, we show 
through water-rock alteration experiments and complementary 
chemical speciation calculations that phosphate-rich fluids are an 
inevitable outcome of alkaline and carbonate-rich ocean water react- 
ing with unaltered carbonaceous chondritic rock. 


Results 


Weconductedacomprehensive survey of 345 Type 3 particles (Methods, 
Extended Data Table 1). Our analysis led to the detection of a rare sub- 
type (nine particles; Methods, Extended Data Table 2) whose spectra 
show highly significant features of phosphate species: in addition to 
the defining mass spectral peaks of Na-salt-rich water ice grains at 23 
and 63 u (ref. 3), there were peaks corresponding to molecular masses 
of 125,165 and 187 u forming a pattern we found to be unique to sodium 
phosphates at high concentrations within these ice grains (Fig. 1 and 
Methods). 

To obtain a quantitative assessment of the composition of these 
phosphate-bearing ice grains, an analogue experiment was used to 
simulate the impact ionization of ice grains viaa laser pulse targeting a 
micrometre-sized water beam”. This technique, laser-induced desorp- 
tion of ions and ionic aggregates”, has previously reproduced accurate 
CDA ice grain spectra*’. Despite the higher amplitude of the (Na;PO,) 
Na’ peak compared with that of (Na,HPO,)Na’, the peak pattern can 
best be reproduced by dissolving Na,PO, and Na,HPO, at an inverse salt 
molar mixing ratio between the two salts of 1:2.5-1:25 (Fig. 2), indicating 
the dominance of monohydrogen phosphate (Na,HPO,) over phosphate 
(Na,PO,) in these ice grains. The dominance of Na,HPO, is substanti- 
ated by the presence of (NaPO;)Na’* (125 u), a charged molecular meta 
phosphate aggregate (Methods, Extended Data Fig. 4) that formed 
from Na,HPO, in our analogue experiments, but not from Na,PO, 
alone (Methods, Extended Data Figs. 1 and 2). With the molar ratio of 
Na,HPO,:Na,PO, kept above 2.5:1.0, a concentration ranging from 0.05 
to 0.60 M of total phosphate salts then gave a close match with CDA 
spectra of phosphate-rich grains (Methods). 

We also tested for the presence of phosphites (for example, Na,HPO,, 
where P hasan oxidation state of +III compared with +V in the orthophos- 
phates Na,PO, and Na,HPO,). Because the acid dissociation constant 
(pK,) of H,PO, is about 6.2 at O °C (ref. 24), this species is unlikely to 
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Fig. 2 | Spectrum from the LILBID analogue experiment reproducing the 
features inthe CDA spectrum. An aqueous solution of 0.420 MNa,HPO, and 
0.038 MNa,PO, was used. All major characteristics of the CDA spectrum of 
phosphate-rich grains (Fig. 1) are reproduced at the higher mass resolution of 
the laboratory mass spectrometer (roughly 700 m/Am). Note: this solution is 
not equivalent to the inferred ocean concentration. To derive the latter 
quantity, the concentration determined in these P-rich grains must be 
averaged over the entire dataset of salt-rich ice grains (see below). 


be relevant to Enceladus’s alkaline ocean?” Raising the concentration 
of Na,HPO, in the phosphate mixture to 5 mM and above led to a clear 
mass line at 149 u due to (Na,HPO,)Na* (Methods and Extended Data 
Fig. 3). However, this peak was not observed in CDA spectra (5.7-5.8 us 
in Fig. 1), thus providing an upper limit on phosphite concentrations 
in these ice grains that is one to two orders of magnitude below the 
required phosphate concentration. 

We now follow the approach established in Postberg et al.?, which 
assumes that the total ensemble of salt-rich grains emitted currently by 
the plume is the best representation of the composition of the source 
water. Thus, to infer the concentration in the ocean we calculate the 
average phosphate concentration in all Type 3 grains. We ‘dilute’ 
the concentration in the nine spectra of the new phosphate-bearing 
subtype (2.6 + 0.9% by number) with the 336 spectra of other Type 3 
grains—which are mostly dominated by NaCl or Na,CO, (ref.3)—and do 
not show phosphates at detectable concentrations ([ZPO,° ] <5 uM). 
Including the +0.9% s.e.m. from counting statistics and the inferred 
concentrations of total dissolved phosphate ([ZPO,° ] = 0.05 - 0.60 M; 
[ZPO] = [PO,?] + [HPO,” ] + [H,PO, ], where [X] denotes the molar 
concentration of species X), this yields a phosphate concentration 
range of 0.8-21 mM for the Enceladean ocean. From our analogue 
experiments, the upper limit inferred from the nondetection of 
phosphites lies at least one order of magnitude below phosphate 
concentrations. This indicates a clear predominance of orthophos- 
phates as the most abundant P-bearing species in the Enceladean 
ocean. 

To interpret the CDA detection of phosphates we performed hydro- 
thermal laboratory experiments with powdered samples from a car- 
bonaceous chondrite (Jbilet Winselwan (CM-2 type))” as an analogue 
of Enceladus’s core material”. The main purpose of the hydrothermal 
experiments was to elucidate the geochemical controls on phosphate 
solubility under conditions relevant to Enceladus. The chondritic 
material was exposed to two different Na-carbonate-enriched solu- 
tions (Run no. 1: initial ZCO, = 0.2 M; Run no. 2: initial SCO, =1.0 M), 
where the total dissolved inorganic carbon concentration is given by 
XCO, = [CO; a4] + [HCO; ] + [CO;?] + [NaCO, ] + [NaHCO,]. The start- 
ing solution also contained 0.5 M XNH; = [NH; a] + [NH,"], a plausi- 
ble concentration for Enceladus’s ocean’. An elevated temperature 
(150 °C) was used in the experiment to accelerate reaction rates. 
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Fig. 3 | Phosphate concentrations derived from hydrothermal experiments. 
Experimental results (squares for Run no. 1 and circles for Run no. 2, see text) 
of P concentrations in hydrothermal fluids at 150 °C and 30 MPa are shown as a 
function of CO, concentration. Numbers attached to data points are calculated 
in situ pH (Methods) and reaction times in hours. Solid and dashed curves 
represent results from thermodynamic equilibrium calculations (Methods) for 
hydroxyapatite-calcite and whitlockite-calcite buffer systems, respectively, 
evaluated at different pH values. 


This allowed us to study key aspects of the water-rock interactions 
occurring in Enceladus at time scales accessible in the laboratory of 
up to about 1,000 h. Fluid samples were collected periodically from 
the reaction cell during the experiments to determine pH and fluid 
composition, including the concentrations of P and XCO, (Methods). 
The initial [£2CO,] adopted for Run no. 2 (1.0 M) was higher than the 
value inferred for the Enceladean ocean’, but a large fraction of initial 
XCO, was removed during the experiments via formation of carbon- 
ate minerals (Methods), leading to more realistic levels of carbonate 
species (Fig. 3). 

Our experiments show that the dissolved phosphate concentration 
increases with in situ measured ZCO, (Fig. 3). This happens because an 
increase in CO,” concentration (or, more correctly, aqueous activity) 
leads to dissolution of Ca phosphate minerals with the formation of sta- 
ble Ca carbonates (that is, calcite)”. In comparison with thermody- 
namic equilibrium predictions, we find that the hydroxyapatite-calcite 
buffer system (hydroxyapatite + 5CO,” + H* + calcite + 3PO,* + H,O) is 
generally applicable to our experimental results (Fig. 3). The relevance 
of this buffer is further supported by chemical analysis of the solid 
residues after Run no. 2, which demonstrates the presence of coexisting 
Ca phosphate with high Ca:P ratios (probably a mixture of hydroxyapa- 
tite and Ca-/Mg-whitlockite) and calcite (Methods, Extended Data 
Figs. 5-7, Extended Data Table 4). One data point in Fig. 3 (CO, = 0.01 M 
in Runno. 1) lies significantly above the hydroxyapatite-calcite buffer, 
appearing to be more consistent with the whitlockite-calcite buffer 
(whitlockite (B-Ca,(PO,),) + 3CO,” < calcite + 2PO,° ). We suspect that 
this change in behaviour is due to unusually high concentrations of Ca” 
and Mg” in Runno.1(0.8 and 0.3 mM, respectively; Methods, Extended 
Data Fig. 8, Extended Data Table 3), which may promote the precipita- 
tion of whitlockite™. This interpretation is supported by an electron 
probe microanalysis of the solid residue (Methods, Extended Data 
Table 4). Alternatively the CM chondrite, Jbilet Winselwan (Extended 
Data Figs. 9 and 10), used in our experiments is known to contain 
impact-induced, thermally metamorphosed materials”, which could 
make some of the results less relevant to Enceladus. By contrast, buffer 
systems involving other phosphates (for example, Mg-whitlockite and 
merrillite) and carbonates (for example, dolomite) cannot explain the 
experimental results (Methods). 


Figure 4a shows the modelled ZPO,” asa function of CO, for water- 
rock interactions adopting the two buffer systems inferred from our 
high-temperature experiments, but now transferred to conditions 
corresponding to the cold seafloor of Enceladus (0.1°C and 8 MPa). 
TheP concentrations determined from the CDA measurements in this 
work (ZPO,° = 0.8-21 mM), the independently derived constraint on 
ICO, (0.04-0.20 M)””8, as well as an alkaline pH (9.0-10.5)*, are all fully 
consistent with these two buffer systems. The hydroxyapatite-calcite 
buffer is applicable for the entire range of ZCO, (Fig. 4a), corresponding 
to PO,” of about 0.8-12 mM. The whitlockite-calcite buffer is appli- 
cable for low ZCO, of about 0.04-0.07 M corresponding to ZPO,” in 
the range of about 7-21 mM at 0.1 °C. 

Figure 4b shows the temperature dependence of ZPO,/ in fluids 
controlled by the hydroxyapatite-calcite buffer. The effects of colder 
fluids, as used in the experiments, down to oceanic temperatures near 
0 °C can be seen. The phosphate concentration increases with lower 
temperature at constant pH because hydroxyapatite (and other reac- 
tants) becomes less stable relative to calcite (and other products) at 
lower temperatures. To achieve XPO,” inthe range 0.8-21 mM, as con- 
strained by the CDA measurements, the temperature of the water-rock 
reactions that provide phosphate in Enceladus would need to be 
80 °C or less if the pH is 9.5, and below 65 °C if the pH is 10.5 (Fig. 4b). 
Higher-temperature fluids could contain sufficient phosphate only 
if the fluid pH were lower. If the whitlockite-calcite buffer controls 
XPO, in Enceladus’s ocean, the temperature could be as high as around 
180 °C (Extended Data Fig. 11) for a maximum ZCO, of 0.2 M. Inthe 
case in which high-temperature interactions with seafloor rocks pro- 
vide the oceanic phosphate, its concentration in the fluids should be 
preserved on cooling to near subzero temperature within the cold 
ocean. This scenario would, however, be relevant only if the time scale 
of cycling the ocean through hydrothermal systems is shorter than that 
of low-temperature dissolution of phosphate minerals at the seafloor. 

Given the consistency of the measured XPO, , XCO, and pH of 
Enceladus’s ocean, we suggest that thermodynamic equilibrium 
between Ca phosphate (for example, hydroxyapatite) and calcite is 
probably achieved in both the cold global ocean and hydrothermal 
fluids of moderate temperature. 


Discussion 


The CDA detection of ice grains with high concentrations of orthophos- 
phates indicates that phosphorus is readily available at the top of 
Enceladus’s ocean (that is, the plume source region). Even with a con- 
servative margin, our estimate indicates concentrations in the order 
of at least hundreds of micromolar, several 100-fold the average phos- 
phate abundance in Earth’s oceans”. The heterogenous distribution of 
salts inthe emitted ice grains might suggest some level of heterogeneity 
in the ocean and that actual concentrations could be locally higher 
or lower than the bulk number calculated in this work. Alternatively, 
yet unknown mechanisms might separate different salts when ocean 
material is transported upwards in the vents. However, even without 
knowing the details of the subsurface processes, we assume that the 
total ensemble of salt-rich grains emitted by the plume is a reasonable 
representation of the average ocean composition near its surface. 

A phosphate:phosphite ratio of at least around 10:1, as implied by 
our results, is consistent with the much greater thermodynamic sta- 
bility of phosphates relative to phosphites under Enceladean ocean 
conditions”. In addition, the nondetection of phosphites suggests that 
aqueous alteration of iron phosphides (the most probable accreted 
form of P) is no longer a significant source of P for Enceladus’s ocean. 
This process would produce metastable phosphites”. The inferred 
lack of contemporary alteration of phosphides is consistent with the 
presence of a low-density rocky core that has already experienced 
extensive water-rock interaction”"****», Phosphites could have been 
produced when accreted phosphides first reacted with liquid water, 
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Fig. 4| Comparison of observed and calculated concentrations of ZPO,” 
in fluids affected by water-rock reactions within Enceladus. a, Relation 
between ZPO,’ and XCO, ata temperature of 0.1°C for the hydroxyapatite- 
calcite buffer system (solid lines) and the whitlockite-calcite buffer system 


(dashed lines). Constraints on ZCO, obtained in previous studies”: are 


but then might have been cycled through hydrothermal systems where 
they would have been converted to more stable phosphates. 

The time scale of our experiments is much shorter than those of 
water-rock reactions in Enceladus, giving rise to uncertainties in the 
exact phosphate/carbonate minerology inside the icy moon. Recent 
work” estimated [ZPO,° ] in Enceladus’s ocean at low temperatures, 
as achieved by interactions involving thermodynamically stable phos- 
phate and carbonate minerals (for example, merrillite and dolomite). 
The upper end of the predicted range from that study also agrees with 
the ocean concentration of phosphates (0.8-21 mM) derived here if 
the ocean pH is 8.0-10.5. This is consistent with the range of 9.0-10.5 
adoptedinthis work. A notable difference is that the range given by ref.17 
also permitted [ZPO,° ] to reach much lower values—specifically at pH 
above 10.0—than those based on the hydroxyapatite-calcite buffer in 
Fig. 4a. The probable reason for this is that, under the strongly carbon- 
ated conditions modelled by ref. 17, hydroxyapatite was more soluble 
and thus merrillite took over the role of the P-controlling phosphate 
mineral. It is important to note that both systems produce results in 
good agreement with the CDA measurements and, in both cases, the 
driver enabling the abundant availability of phosphate is the high 
observed concentration of dissolved carbonate species, which shift 
phosphate-carbonate mineral equilibria toward dissolution of solid 
phosphates into Enceladus’s ocean. 

Anearlier study by Lingam and Loeb” suggested that the availability 
of phosphorus would be the bottleneck of bio-essential elements on 
Enceladus and other icy ocean worlds with potential hydrothermal 
activity and without dry land. Steady-state concentrations “likely lower 
than the corresponding value on Earth by a few orders of magnitude” 
would markedly reduce the prospects for life. Indeed, of the six ele- 
ments—C,H, N, O, P and S—that are generally considered to be critical 
ingredients for life based on water and organic chemistry”, phosphorus 
is cosmochemically the least abundant and has not previously been 
detected at any of the ocean-bearing moons in the Solar System. How- 
ever, the results presented here demonstrate that Enceladus instead 
has a high availability of dissolved P, which is thus extremely unlikely 
to bea limiting factor in the survival of putative life on Enceladus—and 
perhaps also on other ocean worlds that reside beyond the CO, snowline 
inthe Solar System (that is, the distance from the Sun beyond which CO, 
isin a solid (icy) state and is available as a planetary building material). 
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indicated by the blue shaded area. The area highlighted in pink represents the 
range of ZPO,° constrained in this study from CDA data. b, Dependence of 
IPO% ontemperature for the hydroxyapatite-calcite buffer and different 
values of pH and ZCO.. A similar diagram for the whitlockite-calcite buffer can 
be found in Extended Data Fig. 11. 


These wider implications are based on the fact that interactions 
between chondritic rocks and CO,-rich fluids at a water:rock (W:R) 
ratio of around 1.0 would lead to a pH of about 10.0 (refs. 16,36), in 
which dissolved phosphate concentrations tend to maximize” (Fig. 5). 
This is supported by recent findings for asteroid Ryugu, whose parent 
body is thought to have formed beyond the CO, snowline”. Given the 
high solubility of Na-bearing phosphate”, its common occurrence 
in Ryugu samples is qualitatively consistent with our conclusion of 
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Fig. 5| Sensitivity of ZPO,’ of the hydroxyapatite-calcite buffer system 

to pHand XCO, at 0.1°C and 8 MPa. The case of XCO, = 0.2 M (red curve) 
corresponds to the ocean water of Enceladus’ and, perhaps, to icy bodies 
beyond the CO, snowline that originally contained CO, ice. The typical pH of 
such XCO,-rich seawater is around 9.0-11.0 at W:R of approximately 1.0, owing 
to high concentrations of dissolved carbonates'**°. The case of ZCO, = 0.01 M 
(purple curve) might correspond to ocean water of icy bodies within the CO, 
snowline that originally contained no CO, ice. According to ref. 16, ZCO, in fluid 
becomes inthe order of 0.01M by interactions between Cl chondritic rocks and 
pure water at 25 °C, around 1 bar and W:R approximately 1.0. The typical pH of 
such CO,-poor water would be around 11.0-13.0 (refs. 16,36). Red and blue 
shaded areas denote probable pH ranges for waters of icy bodies beyond and 
within the CO, snowline, respectively'***. 


phosphate-rich waters beyond the CO, snowline. Enceladus’s ocean 
could be a harbinger of high phosphorus availability in subsurface 
oceans across most of the outer Solar System. On the other hand, if 
an icy body formed within the CO, snowline both low [ZCO,] and high 
pH (approximately 13.0 at W:R around 1.0) would probably limit 
dissolved phosphate concentrations (Fig. 5). 

The stark contrast between earlier modelling and our results might 
be due to modelling assumptions that were based on scaling fluxes of 
the P cycle on the modern Earth to Enceladus”, not considering the 
fundamental differences between Earth and ocean-bearing moons. 
The most important differences are the much higher concentration 
of carbonate species in alkaline ocean water and the probable pres- 
ence of unrecycled, equilibrated rocks at the seafloor of Enceladus**** 
versus continuous production of more reactive seafloor basalts on 
Earth. Regardless of these theoretical considerations, with the finding 
of phosphates the ocean of Enceladus is now known to satisfy what is 
generally considered to be the strictest requirement of habitability. 
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Methods 


CDA data analysis 

Short description of the CDA chemical analyser subsystem. The 
chemical analyser is the subsystem of the CDA” that provides compo- 
sitional information about an impacting dust particle. Depending on 
the trajectory of the particle, it either hits the central rhodium target 
(chemical analyser target (CAT), diameter 0.17 m), the surrounding 
gold target (diameter 0.41 m) or the inner wall of the instrument. This 
work deals only with impacts onthe CAT. If a dust particle impacts the 
CAT with sufficient energy, it is totally vaporized and partly ionized, 
forming an impact plasma of target and particle ions together with 
electrons, neutral molecules and atoms. The instrument separates this 
plasma, the positive component of whichis linearly accelerated towards 
a multiplier about 19 cm away, used to generate a ToF spectrum. The 
spectrometer is sensitive to positive ions only. The mass resolution 
(m/Am), derived from laboratory experiments with the instrument, 
depends onthe atomic masses of the ions. At 1 u, m/Amis 10, increasing 
to30 m/Amat100 u and upto 50 m/Amat 190 u, although these values 
vary markedly with impact conditions. 

Acquisition ofa spectrum can be triggered by charge thresholds being 
exceeded on either the target or the multiplier ion detector (whena 
certain abundant ion species arrives there, typically H*, H,0* or Na*). In 
the case ofa multiplier-trigger, the spectra show only ions with masses 
higher than the triggering ion species; the triggering species and those 
with lower masses do not appear inthe spectrum. For this work, all P-rich 
spectra were triggered by impact. The spectra are logarithmically ampli- 
fied, digitized at eight-bit resolution and sampled at 100 MHz for a period 
of 6.4 ps after the trigger. The recording period of the high-rate sampling 
modeallows the detection ofions with mass of upto 185-200 u, assuming 
that instrument recording is triggered by the impact itself. 

Because ToF is proportional to the square root of the mass:charge 
ratio of ions, its spectrum in an ideal case also represents a mass spec- 
trum for identical ion charges. The ions created by impact ionization 
in the impact speed regime considered here (roughly 10 km s”) are 
almost exclusively singly charged. Unfortunately, ToF is also influ- 
enced bythe broad distribution of initial ion velocities, slightly varying 
flight paths and plasma-shielding effects*°. For that reason, species of 
identical mass are distributed over arange of sampling points andthe 
mass resolution drops below integer values, usually at 20-30 u and 
above. Therefore, we prefer to show original CDA spectra with anx axis 
showing ToF rather than mass; the latter could mislead the reader to 
intuitively assume an unrealistically high accuracy. 

For a detailed description of the instrument see ref. 10; the calibra- 
tion routine is described in refs. 6,20. For the dataset in this work the 
mass calibration was generally done using sodium (Na*) and sodium 
hydroxide ((NaOH)Na’*) mass lines as reference and a stretch factor 
of a=473 ns. 


CDA dataset. All data used for this work are listed in Extended Data 
Table Land, like all CDA data, are archived on the Small Bodies Node of the 
Planetary Data System (PDS-SBN), at https://sbn.psi.edu/pds/resource/ 
cocda.html. We used 15 time periods between 2004 and 2008 covering 
times of CDA pointing that were favourable for the detection of E-ring 
dust. Times this early in the Cassini tour were chosen because CDA con- 
tamination from salts deposited onthe CDA impact target during deep 
Enceladus plume crossings was then negligible”. Within these periods 
atotal of 7,353 spectra from E-ring grains were recorded. A Lee filter was 
applied to the spectra to improve signal to noise. Of these spectra, 962 
have been categorized as Type 3 (Extended Data Table 1). 

From this set of Type 3 spectra, nine stand out as belonging toa pre- 
viously unknown subtype with a unique pattern of three peaks of m/z 
above 120 u; these are the P-rich spectra discussed above. Individual 
events are listed in Extended Data Table 2. Remarkably, all these spectra 
were triggered on the instrument by the impact charge of the particle 


whereas only 36% of all Type 3 spectra in the dataset were triggered in 
this way (Extended Data Table 1). The remaining 64% of Type 3 spec- 
tra were triggered by a particular ion species reaching the multiplier 
(Methods). By applying binomial statistics, we get a probability that all 
P-rich spectra were impact charge triggered by chance of only 0.01%. 
Therefore, we assume a systematic effect in spectral recognition for 
P-rich grains. Indeed, it is plausible that the spectrum of a P-rich grain 
is not identified as Type 3 if it is triggered by an ion species at the 
multiplier. Type 3 spectra are in most cases triggered by Na’, whichis 
then missing in the spectrum. However, without the Na’ line and inthe 
absence of other typical Type 3 features suchas (NaCl),,Na* or (Na,CO;) 
Na‘, P-rich spectra are not identifiable as being Type 3 following the 
criteria normally used, which focus on spectra dominated by chloride, 
carbonate and hydroxide peaks*”. 

We therefore used only those 345 spectra of the dataset that—like 
the new phosphate-rich spectra—have been triggered by impact 
(Extended Data Table 1) as a reference to calculate the frequency of 
the new phosphate-bearing subtype within the family of Type 3 grains 
and, with that, the average phosphate concentration in salt-rich Type 3 
grains. Itis noteworthy though, that even in the unlikely case that there 
was no such systematic selection effect and we use the entire dataset for 
the inference of average phosphate concentration, the lower limit will 
still be ZPO, above 0.3 mM, more than 100-fold the mean concentra- 
tionin Earth’s oceans, and the conclusions of this work are unaffected. 


Laboratory methods simulating CDA mass spectra 
Identification of sodium phosphates. The largest peaks in all nine 
CDA spectra always appear at 23 and 63 u, indicative of Na* and (NaOH) 
Na’ cations that are the defining features in all salt-rich spectra from 
grains emitted by Enceladus (Type 3; ref. 3), thus these spectra contain 
sodium salts in an alkaline matrix. We used the other three characteristic 
peaks, at 125, 165 and 187 u (allowing a divergence of +1 u that reflects the 
mass uncertainty of CDA”), to perform a rigorous database search— 
in both our own database base“ as well as several external databases 
(for example, NIST Chemistry WebBook, MassBank North America and 
MassBank North Europe (https://massbank.eu/MassBank)). We found 
phosphates to be the only possible class of Na-rich compounds that 
produce a pattern that isin even remote agreement with the observed 
cationic peak pattern. After that we widened the accepted mass range 
to +2 u away from the observed masses and relaxed the restriction 
to Na-rich compounds to determine whether there are other species 
(without sodium) forming a characteristic mass line pattern similar to 
the three mass lines (125, 165 and 187 u). The only match remotely close 
was methyl-2,2-dimethyl-3-(2,2-dichloroethenyl) cyclopropanoate, 
showing prominent peaks at m/z 61(-2), 127(+2), 163(-2) and 187. 
The mass spectrum can be found at ref. 42. 

It is questionable whether a sufficient concentration of this type 
of chlorinated organic ester would be stable in Enceladus’s alkaline 
ocean. Even if we assume that this unlikely substance could be mixed 
with sodium salts in an ice grain, the spectrum has additional features 
not matching the CDA spectrum—most prominently a peak at m/z 91, 
which is absent in all nine CDA spectra. Similarly, any combination of 
exotic organic compounds we found that remotely reproduces the five 
prominent CDA peaks does have anumber of additional peaks that are 
absent in CDA spectra. 

In conclusion, sodium phosphates are the simplest and—to the best 
of our knowledge—the only explanation for all the peaks observed in the 
CDA spectra (Figs. land 2). These salts have already been predicted to be 
relatively abundant in Enceladus’s ocean based on recent geochemical 
modelling”, and the observed peak pattern can be readily reproduced 
in mass spectra from our analogue experiment, as described in the 
following section. 


LILBID experiments simulating CDA spectra. We used the laser- 
induced liquid beam ion desorption (LILBID) facility at Freie Universitat 


Berlin to simulate CDA spectra. The experimental setup is described 
in detail in ref. 22 and here we provide only a brief overview. A 
micrometre-sized water beam is irradiated by a pulsed infrared laser 
(20 Hz, 7 ns pulse length) at a wavelength of 2,840 nm and variable laser 
energy. On absorbing the laser energy, the water beam explosively dis- 
perses into neutral and charged macroscopic, molecular and elemental 
fragments. Cations or anions, depending on the instrument polarity, 
are accelerated and analysed in a ToF mass spectrometer. Detected 
signals are amplified, digitized and recorded with aLabVIEW-controlled 
computer. Typically 300-500 individual spectra are co-added and 
averaged during one effective measurement. The laboratory setup is 
calibrated before every measurement using a 10° M NaCl solution”. 
The recorded mass spectra showed a mass resolution of 600-800 m/Am 
(full-width at half-maximum). 

As discussed in Methods, we found phosphates to be the only Na-rich 
compounds that are in even remote agreement with the observed cati- 
onic peak pattern. Because the appearance of certain peaks and their 
amplitudes are very sensitive to the phosphate species used and their 
concentrations, we individually tested a variety of combinations of 
aqueous solutions of Na,PO,, Na, HPO,, NaH,PO, and Na,HPO,, as well 
as mixtures of these compounds at different concentrations with the 
goal of reproducing the CDA peak pattern. Concentrations were varied 
between 0.01 M and about 5 M for NaH,PO,, 0.5 M for Na, HPO, and 
1.5 M for Na,PO,. 

We found that a mixture of Na,PO, and Na,HPO, most closely repro- 
duced the CDA peak pattern of P-rich spectra (Fig. 1), whereas spectra 
of the pure substances always exhibited differences (Extended Data 
Figs. land 2). 

We also tested for the presence of phosphites (using Na,HPO,, 
in which P has an oxidation state of +III) and hypophosphite (using 
NaH,PO,, in which P has an oxidation state of +I). Both substances cre- 
ated mass lines not in agreement with CDA observations: (Na,HPO;) 
Na’ (149 u) for phosphite and (NaH,PO,)Na* (111 u) for hypophosphite. 
For phosphite we inferred an upper limit concentration of 5 mMtobe 
in agreement with the absence of a mass line at 149 uin CDA spectra 
(Extended Data Fig. 3). 

The phosphate salts used in our experiments were: 

e Na,PO,: Sigma-Aldrich, 96% purity 

¢ Na,HPO,: Roth, over 99% purity 

¢ NaH,PO,°2H,0: AppliChem, over 98% purity 

¢ NaH,PO,°5H,O: Sigma-Aldrich, over 98% purity 
¢ NaH,PO,°H,O: Sigma-Aldrich, over 99% purity 


Theoretical considerations regarding the observed characteristic 
pattern of phosphate ions. We observed the peak sequence at masses 
187, 165 and 125 u as a characteristic pattern indicating the presence 
of phosphorus (P) salts in Enceladus’s ice grains, present in all nine of 
the individual CDA ice grain spectra. Phosphate-containing salts from 
the solvation and extraction of minerals at a basic pH value produced a 
number of anion species in solution that are coupled through a multi- 
step acid-base thermodynamic equilibrium and proton exchange—for 
example, PO, — HPO,” e H,PO, . For the given high salt concentra- 
tions, abundance of counter cations are close to anion abundances in 
solution (Na,PO,°™ ). After ice grain impact, ionic aggregates showed 
additional clustering with abundant cations such as Na‘ resulting in 
characteristic mass peaks due to species including (Na,HPO,)Na* and 
(Na,PO,)Na*. Owing to the basic pH of the solution forming the ice 
grains (or the liquid desorption matrix of the laboratory experiment), 
these two species are expected. 

However, we did not detect a (NaH,PO,)Na* adduct in CDA spectra 
but rather a (NaPO,)Na* peak at 125 u, a charged molecular aggregate of 
ametaphosphate with phosphorus in the same oxidation state (+V) as 
the other orthophosphates. We attribute the absence of (NaH,PO,)Na" 
and the appearance of (NaPO,)Na’ instead to an elimination reaction 
during desorption (LILBID) or particle impact (CDA): 


NaH,PO, + heat/energy > NaPO, + H,O, followed by the addition 
of Na’. 

The fact that the (NaPO,)Na* peak at 125 u is formed from hydrogen 
phosphates (linked via an acid-base equilibrium) in our analogue exper- 
iments, but not from Na,PO, alone, supports this picture. However, it 
should be noted that we could produce a (NaH,PO,)Na* peak (143 u) 
if high concentrations (over around 0.05 M) of NaH,PO, or Na, HPO, 
(Extended Data Fig. 2) were used in the laboratory experiments or if 
the pH value of the solution was forced below about 9.0. We then find 
aclear correlation of the amplitudes of dihydrogen phosphate (143 u) 
with metaphosphate (125 u). In the CDA data the 143 u peak appears to 
be below the noise level of the detector. In our laboratory spectra the 
125 umetaphosphate peak is generally three- to fivefold more intense 
and itis therefore likely that small amounts of (NaH,PO,)Na‘ (143 u) were 
buried in the noise of the CDA spectrum. Such a signature was indeed 
tentatively observed in one of the nine spectra (event 3 in Extended 
Data Table 1). 

To test these ideas we calculated the thermochemical properties of 
the water elimination of NaH,PO, + heat/energy > NaPO, + H,O with 
the Jaguar quantum chemistry programme package at the PBEO-D3/ 
6-311+G(d,p)/PBE level of theory. Two variants have been calculated: 
(1) with Na* as a counter ion and (2) without Na+ (involving a H,PO, 
anion). The reaction is slightly endergonic in the gas phase and even 
more favoured in/with water. The barrier of the concerted elimination 
reaction appears to be low if water molecules are involved. One may 
even consider it to be a water-assisted reaction, in particular for the 
case of Na* as counter ion (Extended Data Fig. 4). 


Hydrothermal experiments and calculations of chemical 
equilibria 
Experimental methods and procedures. Hydrothermal experiments 
were performed with a Dickson-type hydrothermal apparatus (Toyo 
Koatsu Co. Ltd) used in a previous study®. In this apparatus, water-rock 
reactions occur ina flexible reaction cell that comprises a gold bag with 
atitanium head. The flexible reaction cell was connected to a stainless 
steel sampling tube with a gold-coated inner wall, to avoid catalytic 
reactions. The reaction cell was heated at 500 °C for 3 hin air before 
each experiment to remove potential contamination by organic matter. 
The reaction cell was placed in a stainless steel autoclave, which was 
held at a pressure and temperature of 30 MPa and 150 °C, respectively, 
during the experiments. The reaction cell was pressurized by the addi- 
tion of surrounding water into the autoclave using a hand pump. The 
autoclave was heated with an external heater. Due to the flexibility of 
the gold bag, the reaction cell deforms plastically without fracturingin 
response to the addition of surrounding water. This allowed us to under- 
take online sampling of fluid without significant changes in pressure. 
The starting rock material was prepared by powdering a few pieces of 
the commercially available Jbilet Winselwan carbonaceous chondrite 
(CM-2 type)” in an agate mortar. The starting solutions were prepared 
by dissolving NaHCO; powder and NH; solution in ultrapure water such 
that the initial ZCO, was 0.2 M for Run no. 1 and 1.0 M for Run no. 2, 
and XNH, was 0.5 M. One way to achieve high XCO, in an Enceladus 
experiment is to introduce gaseous CO, in the reaction cell and dis- 
solve it by pressurization. However, this is technically challenging in 
our experimental system and we cannot control the [ZCO,] and pH of 
the starting fluids. Therefore, we achieved high XCO, by the addition 
of NaHCO, powder. We added 1.0 M Nain Run no. 2, higher than Na 
concentrations in Enceladus’s ocean’. If Na phosphate (for example, 
Na struvite) had precipitated in such Na-rich fluids of Run no. 2, this 
would have affected the buffer system that controls ZPO, in fluids. 
Although the temperature and pressure dependences of the thermo- 
dynamic constants of Na struvite have not been experimentally deter- 
mined, the saturation index (log([Na*][Mg”*][HPO,” ]/Kya-struvite) Where 
Kya-struvire denotes the equilibrium constant from ref. 17 of Na struvite 
at 1 bar and 25 °C, is calculated to be below -1.6 using the measured 
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Na, Mg and ZPO, concentrations and in situ pH. Because the fluid Mg 
concentrations of Run no. 2 were below the detection limit (Extended 
Data Fig. 8b), Na struvite was likely to have been undersaturated in 
the experiment. The presence of Ca sulfate in the starting rock materi- 
als (CM chondrite) also does not affect our conclusion regarding the 
hydroxyapatite-calcite buffer system. If the hydroxyapatite-gypsum 
(or hydroxyapatite-anhydrite) buffer controls ZPO,”, [ZPO,° ] 
should dramatically increase with [SO,” ]; however, this did not occur 
(Extended Data Table 3). The initial amounts of starting rock powder 
and solution were around 7 g and 10 ml, respectively. In each experi- 
ment we collected several aliquots (about 1.2 ml) of fluid for analysis. 
Assuch, water:rock mass ratios in the reaction cell changed from around 
1.4-0.9 due to the removal of multiple fluid samples. 

The collected fluid samples were analysed by inductively coupled 
plasma atomic emission spectroscopy (ICP-AES: SPS5510, Hitachi 
High-Tech) for dissolved Na, K, Mg, Ca, Siand Al content; by ICP mass 
spectrometry (ICP-MS: Agilent 8000) for dissolved Fe content; by 
ion chromatography (ICS-1600, DIONEX) for dissolved Cl and SO, 
content; and by gas chromatography (GC-2010, Shimadzu) for £CO, 
content (Extended Data Table 3). PO,? content was analysed by 
both tandem ICP-MS (ICP-MS/MS: Agilent 8000, Agilent) and ion 
chromatography. The former methodology quantified dissolved 
total P in fluids at high sensitivity, including phosphate and poten- 
tially phosphite, whereas the latter directly measured XPO,° with 
moderate sensitivity. The results of the two analyses showed that, 
basically, all dissolved total P in the fluids was phosphate (Extended 
Data Table 3). As such, we used ICP-MS/MS results to quantify ZPO, 
inthe fluid samples. In situ pH values at 150 °C during the experiments 
were calculated using the geochemical code PHREEQC v.3 (ref. 44), 
based on the measured pH at 25 °C and dissolved elements, molecules 
(phosphate and sulfate) and CO, in the diluted fluid samples. To 
calculate the in situ pH of the alkaline, Na-carbonate-rich fluids, we 
constrained the charge balance from the measured pH at 25 °C in two 
ways. First, Na was used to compensate for charge imbalance due to 
analytical uncertainties. Second, XCO, was used to compensate for 
charge imbalance. We adopted the median value between the two 
values as the in situ pH, and the difference between the median and 
upper/lower values as the error. We also conducted mineralogical 
and chemical analyses of solid residues after the experiments using 
an X-ray diffraction spectrometer (XRD: MiniFlex600, Rigaku) anda 
scanning electron microscope with an electron probe microanalyser 
(EPMA: JXA-8530F, JEOL). 

The XPO," concentrations controlled by the hydroxyapatite-calcite 
and whitlockite-calcite buffer systems were calculated with the equilib- 
rium constants obtained from the SUPCRT92 programme”. In the cal- 
culations, the thermodynamic properties of whitlockite (B-Ca,(PO,),) 
and hydroxyapatite were taken from Robie et al.*°. 


Elemental composition of carbonate and phosphate. Extended Data 
Figs. 5 and 6 showtypical secondary electron microprobe images and 
elemental mapping for thick sections of solid residues after the experi- 
ments (Run nos. land 2). Inthe solid residues we found large grains of 
precipitated Ca carbonate (CaCO,, size over 100 um; Extended Data 
Fig. 5 and Extended Data Table 4). Because the typical size of Ca carbon- 
ate in the sample before the experiments was 10 um or less (Extended 
Data Fig. 7), the large Ca carbonate grains found after the experiments 
are highly likely to have been precipitated during the experiments 
through XCO, sequestration. This view is supported by decreases in 
XxCO, in fluids (Extended Data Fig. 8). Our XRD analysis shows that the 
precipitated Ca carbonates are most likely calcite in both of the experi- 
ments (Run nos. 1 and 2; Extended Data Fig. 10). 

Toinvestigate the composition of the phyllosilicate phase before and 
after the experiments, we performed EPMA analysis for the phyllosili- 
cate matrix around Ca phosphate grains. Extended Data Fig. 9 shows 
aternary diagram with the results of EPMA analysis for phyllosilicate, 


Ca phosphate and Fe oxides inthe samples before and after the experi- 
ments. Phyllosilicates in the starting materials were mostly serpentine; 
however, mixtures of serpentine and saponite appeared in the samples 
after the experiments. Extended Data Fig. 9 indicates a linear trend ofa 
mixing line of the compositions between Ca phosphate and serpentine 
(and saponite mixtures). 

Extended Data Table 4 presents the results of spot analysis for Ca 
phosphate and carbonate found in the solid residues after the experi- 
ments of Runnos. land 2. At the micrometre scale, Ca phosphate min- 
erals usually coexist with other secondary minerals such as serpentine 
and saponite (Extended Data Fig. 6), and thus the results of the spot 
analysis might not be explained by the simple chemical composition 
of analteration mineral. As reported in Extended Data Table 4, Grain C 
with the highest P content had high total mass fractions of around 
95-100 wt%, suggesting that the contribution of calcite to Ca:P atomic 
ratios is small (less than 2-3% of Ca mole abundance), whereas serpen- 
tine and saponite would coexist with Ca phosphate of Grain A, Vein B 
and Grain D. These hydrated secondary minerals can reduce the total 
mass fractions of spot analyses by 5-10% owing to the presence of H,O 
and OH. Because the total mass fractions for Grain A, Vein Band Grain D 
are also high in general (over 87 wt%; Extended Data Table 4), calcite 
is presumed to be rare in Ca phosphate grains, suggesting that the 
contribution of calcite to Ca:P ratios is small. 

Grain C of Runno. 2 has very high P:Si ratios (Extended Data Table 4c). 
The measured Ca:P mole ratios are 1.32-1.37. Grain D of Run no. 2 has 
relatively high P:Si ratios, above 1.0 (Extended Data Table 4d). The 
Ca:P mole ratios of Grain D range from 1.29 to 1.34. The Ca:P ratios 
of hydroxyapatite, whitlockite [B-Ca,(PO,),] and Mg-whitlockite 
[Ca,MgH(PO,),] are 1.67, 1.50 and 1.29, respectively. Ca phosphates 
collected after Run no. 2 (Grains C and D) are inferred to consist of 
mixtures of multiple Ca phosphates. Due to the trace amounts of Ca 
phosphates in the solid residues, these could not be identified in the 
XRD spectra (Extended Data Fig. 10) 

In Ca phosphate grains of Run no. 1 (Grain A and Vein B), P:Si mole 
ratios are generally low (roughly 0.1-1.0) (Extended Data Table 4a,b), 
suggesting that Ca phosphates coexist with serpentine and saponite. 
There is a relatively wide variation in the Ca:P mole ratios of Grain A and 
Vein B, from approximately 1.1to 1.5 (Extended Data Table 4a,b). A high 
Ca:P ratio of around 1.5 requires the presence of hydroxyapatite and/or 
whitlockite. To explain the low Ca:P ratios (below 1.29) of the analysed 
spots on Grain A and Vein B, the presence of Ca-poor phosphate, such 
as NH,MgPO,, would be needed. 

Our results showing high C:P ratios (1.32-1.37) for Run no. 2 and low 
Ca:P ratios (1.25-1.27) for Run no. 1 are consistent with our interpreta- 
tion of the buffer system that controls ZPO,” concentrations in the 
fluid samples. The Ca phosphate grains (C and D) of Run no. 2 contain 
high fractions of phosphate with high Ca:P, such as hydroxyapatite. 
This observation is consistent with the control of ZPO,” concentra- 
tions in fluid samples of Run no. 2 by the hydroxyapatite-calcite buffer 
(Fig. 3). In contrast, the observation of Ca phosphate grains with low 
Ca:P for Run no. 1 (for fluid with ZCO, = 0.01 M and reaction time 
1,000 h) is consistent with the control of ZPO,” concentrations by 
Ca phosphate with lower Ca:P ratios, such as whitlockite, rather than 
control by hydroxyapatite (Fig. 3). This fluid sample was collected in 
the final sampling of Run no. 1 (Extended Data Fig. 8a), and the solid 
residue was collected after this final sampling. We consider that Ca 
phosphate with relatively low Ca:P, such as whitlockite [B-Ca,(PO,),], 
Mg-whitlockite and NaMgPO,, would have been precipitated before 
the final sampling in response to an increase in Mg concentration in 
the fluid (Extended Data Fig. 8), and that the whitlockite-calcite buffer 
would have controlled the P concentration of the fluid at the end of 
experimental Run no. 1. 


Chemical equilibrium considerations regarding phosphate buff- 
ering systems. Extended Data Fig. 11 shows the calculated results of 


[ZPO] in fluids as a function of temperature for different pH and 
[ZCO,] values, assuming the whitlockite-calcite buffer system (see 
above for the equation). Although our experimental data suggest that 
the whitlockite-calcite buffer could have controlled ZPO,” concentra- 
tions inthe fluid at the end of Run no. 1, we consider that this could have 
been an experimental artefact reflecting transient states of supersatu- 
ration of Mg-bearing phases. In the experiments, the dissolution of 
anhydrous primordial minerals (for example, olivine and pyroxene) 
would have provided Mg and Ca for the fluids. These Mg and Ca ions 
would have been removed by reactions with CO,” , forming carbon- 
ate minerals (Extended Data Fig. 5). However, at the end of Run no. 1, 
XCO, in the fluids was depleted (Extended Data Fig. 8) owing to CO, 
sequestration via carbonate mineral formation; therefore, dissolved 
Mg and Ca concentrations increased in the fluids (Extended Data Fig. 8). 
In longer-period reactions, excess Mg and Ca ions would have been 
consumed through the formation of serpentine and saponite. Thermo- 
dynamic equilibrium calculations show that Mg and Caconcentrations 
in fluids that interact with carbonaceous chondritic rocks are typically 
low (0.001-0.010 mM each) in alkaline fluids (pH over 9.0) in the pres- 
ence of serpentine and saponite**””. In Enceladus, water-rock reactions 
would be driven closer to acid-base equilibrium because the reactions 
can occur over much longer time scales than apply in our experiments. 
If the ocean of Enceladus were consistent with the expectation of low 
Mg and Ca concentrations (for example, under 0.1 mM each), then 
XPO,” in fluids in contact with chondritic rocks would be controlled 
by the hydroxyapatite-calcite buffer system. 

In this study, we considered hydroxyapatite and whitlockite 
[B-Ca,(PO,),] as phosphate minerals of buffer systems that control 
[ZPO,° ]in fluids. However, previous work has suggested the possi- 
bility of the occurrence of other phosphate minerals in water—-rock 
reactions in Enceladus”, including Mg-whitlockite and merrillite. 
Using the equilibrium constants of dissolution of Mg-whitlockite and 
merrillite at 1 bar and 25 °C (ref. 48), and assuming [Mg”*] as roughly 
10° M in the fluid samples (Extended Data Table 3), PO,°> values 
controlled by the Mg-whitlockite-calcite and merrillite-calcite buffer 
systems (that is, 6 x 10° and 7 x 10°°M, respectively, for ZCO, = 0.1 
M) are about 1/30 to 1/250 that of the hydroxyapatite-calcite buffer 
(that is, 2 x 10° M for ZCO, = 0.1 M). Although the temperature and 
pressure dependences of the thermodynamic equilibrium constants 
of both Mg-whitlockite and merrillite are unknown, the observed 
XPO, values in our experiments are unlikely to be explained by these 
buffer systems. 

We also assumed that calcite is the carbonate mineral phase for 
the buffer system, but dolomite is usually more abundant than cal- 
cite in CI chondrites*®. However, ZPO,” values controlled by the 
hydroxyapatite-dolomite buffer are orders of magnitude higher than 
those controlled by the hydroxyapatite-calcite buffer for XCO, = 0.1M 
and the typical Mg concentration in our experiments ([Mg”*]=10 M~’). 
Therefore, we conclude also that the hydroxyapatite-dolomite buffer 
may not explain the observed PO,” in our experiments. 
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Extended Data Fig. 1| LILBID spectrum ofa 0.1M solution of Na;PO,. The 
(NaPO,)Na* peak (125 u) cannot be reproduced with sufficient amplitude and 
the ratio of (Na,PO,)Na* (187 u) to (Na,HPO,)Na* (165 u) is too high compared 


Mass [u] 


with CDA spectra. Both these disagreements disappear if sufficient amounts of 
Na,HPO, are added to the solution. 
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Extended Data Fig. 2|LILBID spectrum ofa0.1MsolutionofNa,HPO,.Amass notpresentin a mixture of Na,HPO, and Na,PO, (Fig. 2). This motivates our 
line at 143 u ((NaH,PO,)Na’) appears, whichis notapparentinCDAspectraexcept conclusion on a mixture of the two orthophosphates in the ice grains from 
tentatively in spectrum #3 of Extended Data Table 2. This signature at 143 u is Enceladus analysed by CDA. 
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Extended Data Fig. 3 | LILBID spectrum with a 0.5 M phosphate concentrations below 0.005 M. Note, the CDA detection limit of phosphite in 
concentration (similar to Fig. 2) and 0.01 M phosphite (Na,HPO,). The pure water (without phosphates) is approx. 1nM, more than 3 orders of 
phosphite peak Na,HPO,Na’ (149 u) is not in agreement with CDA spectra and magnitude lower. 

only becomes undetectable in phosphate-bearing CDA spectra at phosphite 


(a)-Reaction NaH,PO, --> NaPO; + H20 
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Extended Data Fig. 4 | Jaguar quantum chemistry program package at the 
PBEO-D3/6-311+G(d,p)/PBE level of theory of the elimination reaction from 
dihydrogen phosphate leading to the observed meta phosphate species. 
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Calculated at PBEO-D3/6-311+G(d,p)/PBE level of theory. 


The presence of water significantly reduces the meta phosphate formation 
barrier, specifically in the case with Na‘. 
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Extended Data Fig. 5| Typical Backscattered Electrons (BSE) SEM images were absent from the sample before our experiments (Extended Data Fig. 7). 
and elemental mapping for solid residues collected after experiment XRD analysis (Extended Data Fig. 10) shows that the precipitated Ca carbonates 
Run #1. The presence of large Ca carbonate grains (>100 pm) indicates the are most likely calcite. 


formation of calcite during the experiment. Such large Ca carbonate grains 
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Extended Data Fig. 6| (a) Typical Backscattered Electrons (BSE) SEM 
images and elemental mapping results for Ca phosphate minerals inthe 
solid residues collected after experiment Run #1. Phosphate-rich grains 
witha size of several micrometres and veins with a width of several 
micrometres and length of -100 pm are seen near a large serpentine grain. In 
addition to Caand P, Nais enriched in Ca phosphate minerals. Close-up images 
of phosphate minerals (the bottom two images) indicate that Ca phosphate 
coexists with other minerals (e.g., serpentine and saponite). (b) Same as 


P ` 
EA 
Ca phosphate (Grain C) 


Ca sea (Grain D) 


N 


. ` Ca carbonate 


é l 


Extended Data Fig. 6a, but for experiment Run #2. Phosphate-rich grains with a 
size of several micrometres are observed (see the bottom two SEM images for 
close-up views of phosphate minerals). Close-up images of Grain D indicate that 
Ca phosphate minerals coexist with other minerals, whereas the surface of 
Grain C exhibits relatively uniform composition. The chemical compositions of 
Ca phosphate (Grain A and Vein B) revealed by microprobe spot analysis are 
summarized in Extended Data Table 4. XRD analysis (Fig. 10) shows that the 
precipitated Ca carbonates are most likely calcite. 
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Extended Data Fig. 7 | Typical SEM images and elemental mapping results of carbonaceous chondrite samples before the experiments. The size of Ca 
carbonate grains is typically -10 pm or less. 


Concentration (M) 


0 200 400 


600 800 


1000 


1200 


Time (hours) 


Concentration (M) 


0 200 


400 


600 800 1000 1200 


Time (hours) 


Extended Data Fig. 8 | Temporal variations in the concentrations of 
selected elements (Na, K, Si, Mg, Ca, and P) and dissolved total inorganic 
carbon (ZCO,). Results for (a) Run #1 and (b) Run #2 are shown. The arrows for 
Mg concentrations at 650 and 1000 h of reaction time for Run #2 indicate that 
the concentrations were less than the detection limit (0.036 mM). The 
measured data and calculated in situ pH of fluids are summarized in Extended 


Data Table 3. The decreases in CO, in both experiments over time were caused 
by the formation of carbonate minerals (i.e., calcite), as confirmed by 
microprobe and the results of XRD analysis shown in Extended Data Figs. 5 

and 10, respectively. Where no error bars are shown, the analytical errors are 
smaller than the data points. All measured data, analytical errors, and 
calculated in situ pH of fluids are summarized in Extended Data Table 3. 
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Extended Data Fig. 9 | Ternary diagram of the results of EPMA analysis of the chemical compositions of phyllosilicate and Ca phosphate-rich grains inthe 
samples before and after the experiments, plotted as molar fractions. Green stars indicate the compositions of pure minerals found in the samples. 
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Extended Data Fig. 10 | Comparison of X-ray diffraction (XRD) spectra of intensities are offset for clarity. The primary XRD peaks of minerals that 
the solid samples before the experiments (red), after experiment Run #1 correspond tothe XRD peaks are also indicated. 
(purple), and after experiment Run #2 (blue). a.u., arbitrary units. The XRD 


Article 


10° 


in-situ pH 
— 95 


=PO,°- (mM) 


40° 
0 50 100 150 200 


Temperature (°C) 


Extended Data Fig. 11 | Results of calculated ZPO," in fluids as a function of temperature for different pH and ZCO, assuming the whitlockite-calcite 
buffer system. Pink shading represents the observationally-derived range of PO," from the present study. 


Extended Data Table 1| CDA set of Type 3 spectra used for this work 


Period # | Period from/to | Period from/to | Total Number | Type 3 spectra 
(SCLK) of Type 3 triggered by 


spectra impact 


1 2004-10-27 23:50:16 
2004-10-28 23:25:50 


2005-03-08 20:03:02 
2005-03-10 02:41:35 


2005-06-26 09:35:30 
2005-06-26 18:40:44 


2005-11-26 08:29:18 
2005-11-26 23:31:28 


2005-12-24 04:19:14 
2005-12-25 06:50:35 


2006-01-16 06:21:33 
2006-01-18 14:44:13 


2006-02-25 00:58:00 
2006-02-26 02:10:00 


2006-03-21 10:14:22 
2006-03-22 18:02:44 


2006-12-03 00:24:12 
2006-12-03 06:33:08 


2006-12-15 04:41:46 
2006-12-15 16:59:56 


2007-05-10 19:14:05 
2007-05-10 19:59:44 


2008-05-09 22:56:12 
2008-05-10 02:19:47 


2008-11-16 20:59:47 
2008-11-16 21:25:45 
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Extended Data Table 2 | Events of phosphate-rich ice grain recorded by CDA. Saturn distance is given in Saturn radii 
(equatorial radius R,=60268km) 


Saturn Radial | Impact Speed 
Distance (Rs) | Estimate (km/s) 
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Since the impact speed onto CDA cannot be measured directly for these grains, the relative velocity of the spacecraft to a circular Keplerian orbit provides a rough estimate. 
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Extended Data Table 3 | Temporal variations in dissolved species in fluids (in M), with pH at the calculated in-situ pH at 150°C 
and 30 MPa 


Run #1 
Time in situ F 
ha Bi SCO; p PO, Mg Ca Na si K Al Fe ci S04 
332 na SEO on ae 1.4E-05 116-04 286-01 916-04 186-03 7.1£-06 5.3E-06 ay 
(3E-03) (4E-06) (2E-05) (2E-02)  (5.5E-04) (2E-04) (4E-06) (1.3E-06) - 
936.3 10.2  8.6E-02 27E-05  2.8E-05  6.5E-05 616-05 266-01  9.1E-04 4.2E-03 yar og  17E-06 34E-02 5.7E-02 
3 (0.3)  (2E-03) (8E-07) (1.1E-05) (1.4E-05) (1.0E-05) (2E-02) (5.5E-04) (5E-04) í (4E-08) (3E-04)  (4E-04) 
10076 10:5  13E-02 15E-05  2.1E-05  3.1E-04 7,6E-04 2.4E-01 1.2E-03  4.4E-03 eg ge gg  31E-06 3.2E-02 9.0E-02 
6 0.2) (2E-3) (2E-07) (1.5E-05) (7E-05)  (1.2E-04) (2E-02) (7E-04) (5E-04) (1E-08)  (4E-04)  (5E-04) 
Run #2 
Time In situ P 
thes a 5CO, p PO, Mg Ca Na si K Al Fe cl 50, 
536.2 102 4.4601 25E-04 = -2.5E-04 5.6E-05 5.8E-05 1.26100 146-03 5.6E-03 23 ap og 3.9E-05 3.2E-02 5.0E-02 
2 (0.5) (4-02) (9E-06)  (7E-05)  (4E-06)  (2E-06) (2E-01)  (7E-04) (1.2E-03) <> (2E-08) (5E-04) (1E-03) 
10.2  3.6E-01 1.9E-04  2.0E-04 3.1E-05 9.9E-01  1.5E-03 6.6E-03 1.3E-05 3.9E-02 8.7E-02 
642.0 5) (5E-02) (4E-06) (1.0E-04) <3-7E-05 (1E-06) (1.5E-01) (BE-04)  (1.4E-03) <3-3E-05 (8E-08)  (7E-04) (1E-03) 
10.1  3.0E-01 1.8E-04  1.9E-04 2.7E-05 9.0E-01 1.2E-03 6.5E-03 7.0E-06 3.9E-02 1.1E-01 
977.905) (6E-02)  (6E-06) (96-05) < 3-6E-05 (9E-07) (1.4E-01) (6E-04)  (1.3E-03) < 34E-05 (2E-08)  (5E-04)  (9E-04) 


Analytical errors are shown in parentheses. For P concentrations, the values of “P” are the results of total dissolved P in fluids measured using a tandem inductively coupled plasma mass spec- 
trometer (ICP-MS/MS); and those of “PO,” are the results of 2PO,° measured using ion chromatography (see Methods 3a). In cases where element concentrations were too low to be measured 
(i.e., Mg and Al), the upper limits of concentration (detection limits) are given. 
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Extended Data Table 4 | Results (mole abundance in mol/kg) of electron probe microanalyzer (EPMA) analyses for Ca 
phosphates [Grain A and Vein B for Run #1, and Grains C and D for Run #2 (see Extended Data Fig. 6)], and for Ca carbonate 
minerals for Runs #1 and #2 


a. Grain A {Run #1) (number of spot analyses = 4} 


Si Al Fe Mg Ca Na K P S Ni P/Si Ca/P 
A-1 2.66 0.27 2.25 2.33 432 058 0.01 3.27 0.35 0.17 1.23 1.32 
A-2 2.51 0.22 3.18 2.01 2.59 0.49 0.01 2.06 1.19 0.37 0.82 1.26 
A-3 4.00 042 2.75 2.86 2.79 0.46 0.03 2.46 0.04 0.06 0.61 1.13 
A-4 3.84 0.51 2.04 2.65 3.18 049 0.02 2.33 0.10 0.07 0.61 1.37 


Mean Ca/P {a) 1.27 (0.10) 
b. Vein B {Run #1) (number of spot analyses = 6} 


Si Al Fe Mg Ca Na K P S Ni P/Si Ca/P 
B-1 454 065 3.22 3.26 192 034 0.03 162 013 0.17 0.36 1.18 
B-2 3.90 041 2.77 2.99 2.80 037 0.02 2.30 012 0.12 0.59 1.22 
B-3 3.42 037 3.79 2.41 2.18 026 001 146 0.36 0.21 0.43 1.49 
B-4 447 061 531 293 1.33 038 002 109 0.07 0.12 0.25 1.22 
B-5 5.68 080 3.38 349 094 032 004 078 O05 0.10 0.14 1.20 
B-6 3.34 041 3.2 2.55 2.34 033 002 2.02 na 0.47 0.60 1.16 


Mean Ca/P {a) 1.25 (0.12) 
c. Grain C {Run #2) (number of spot analyses = 3) 


Si Al Fe Mg Ca Na K P S Ni P/Si Ca/P 
C-1 0.53 005 100 0.46 7.56 061 0 5.52 0.03 0.03 10.51 1.37 
C-2 0.49 0.06 0.75 O52 788 1.01 0 5.97 0.01 0.02 12.28 1.32 
C-3 0.50 0.06 0.77 0.66 759 1.21 0O 5.64 0.02 0.01 11.28 1.35 
Mean Ca/P 1.34 
d. Grain D {Run #2) (number of spot analyses = 2} 
Si Al Fe Mg Ca Na K P S Ni P/Si Ca/P 
D-1 2.37 041 2.26 1.70 3.89 062 0.01 3.02 037 0.07 1.27 1.29 
D-2 1.82 034 2.37 148 478 097 001 3.58 0.07 0.04 1.97 1.34 
Mean Ca/P 1.31 
e. Ca carbonate grains {Run #1) (number of spot analyses = 3} 
Si Al Fe Mg Ca Na K P S 
E-1 0 0 0.04 <0.01 10.4 0.01 0 0.01 <0.01 
E-2 0 0.01 0.06 0.03 10.4 0.02 <0.01 0.01 0.01 
E-3 0.09 0.02 0.13 0.38 10.1 0.02 0 0 0.05 


f. Ca carbonate grains {Run #2) (number of spot analyses = 3) 


Si Al Fe Mg Ca Na K P S 
F-1 0 0.01 0.08 0.01 11.01 0.02 <0.01 0.02 <0.01 
F-2 0 <0.01 0.15 0.04 10.33 0.02 0 0.03 0.01 
F-3 <0.01 0 0.10 0.01 10.75 0.03 <0.01 0.04 <0.01 


The mass fractions of Fe and S were calculated assuming their forms of FeO and SO . 


